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In this study we consider the nonlinear radial oscillations exciting in LBV-stars with effective 
temperatures 1.5 • 10 4 K < T eff < 3 • 10 4 K, bolometric luminosities 1.2 • 1O 6 L < L < 1.9 • 1O 6 L 
and masses 35.7M < M < 49.1M . Hydrodynamic computations were carried out with initial 
conditions obtained from evolution sequences of population I stars (X = 0.7, Z = 0.02) with 
initial masses in the range 7OM < Mzams — 9OM . All models show instability against 
t— 5 . radial oscillations with amplitude growth time comparable with dynamical time scale of the 

star. Radial oscillations exist in the form of nonlinear running waves propagating from the 
boundary of the compact core to the upper boundary of the hydro dynamical model. The 
velocity amplitude of outer layers is of several hundreds of km/s while the bolometric light 
amplitude is AMboi < 0.2 mag. Stellar oscillations are not driven by the ft-mechanism and are 
due to the instability of the gas with adiabatic exponent close to the critical value I\ = 4/3 
due to the large contribution of radiation in the total pressure. The range of light variation 
periods (6 day < II < 31 day) of hydrodynamical models agrees with periods of microvariability 
observed in LBV-stars. 
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INTRODUCTION 

The most luminous (L ~ 1O 6 L ) stars represent a group of luminous blue variables (LBV) 
which involve such well known objects as r\ Car, P Cyg, S Dor and the Hubble-Sandage variables 
in galaxies M31 and M33. All LBV-stars are massive population I stars at the early helium 
burning stage and their small number is due to the rapid evolutionary movement across the 
HRD to the higher effective temperatures. LBV-stars are thought to precede the Wolf-Rayet 
evolutionary stage (Maeder, 1983; Langer et al., 1994) though it is not excepted that they 
might also be the supernova progenitors (Gal- Yam et al., 2007; Trundle et al., 2008). 

Photometric variability of LBV-stars can be divided into three distinct types depending 
on the characteristic time scale (Lamers, 1987; Leitherer et al., 1992). In variability of the first 
type (giant eruptions) the bolometric luminosity increases by several magnitudes and the mass 
of the ejected material is as high as ~ M (Humphreys & Davidson, 1994). The energy of 
the giant eruption can be so high that the LBV-star becomes the SN-impostor (Smith et al., 
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2009). There are known two giant eruption events in the Galaxy. One of them was P Cyg in 
the early XVII century (de Groot, 1988) and another was rj Car in the middle of XIX century 
(van Genderen, 1984; Frew, 2004), so that the time scale of giant eruptions is thought to be of 
the order of > 10 2 yr. In variability of the second type (S Dor type) with time scale of ~ 10 yr 
the bolometric light seems do not change, whereas the visual light changes by about 2 mag 
(van Genderen, 2001). During the absence of S Dor type variability one can observe the cyclic 
light variations (microvariability) with amplitude of < 0.2 mag and characteristic time ranging 
from one to several dozens of day. 

The origin of the both giant eruptions and S Dor variability remains still unclear and 
only microvariations are interpreted in terms of stellar pulsations (van Genderen, 1989). This 
assumption is supported by the linear theory of stellar pulsation (Glatzel & Kiriakidis, 1993; 
Kiriakidis et al., 1993; Dziembowski & Slawinska, 2005) which predicts that outer layers of 
LBV-stars are unstable against radial oscillations. This fact seems to be of great importance 
because observational estimates of microvarion periods can provide us with independent method 
of mass determination. 

The goal of the present study is to consider nonlinear radial stellar oscillations of 
microvariability in LBV-stars. Hydrodynamic modeling of this phenomenon is too complicated 
because of the large ampliture oscillations exciting at the boundary of the dynamical instability, 
so that only two reports on this item have been published by now (Dorfi et al., 2000; Dorfi & 
Gautschy, 2002). Below we present the results of hydrodynamic calculations of radial oscilla- 
tions in LBV-stars with effective temperatures 1.5 • 10 4 K < T cS < 3 • 10 4 K and luminosities 
1.2 • 1O 6 L < L < 1.9 ■ 10 6 L Q . Initial conditions for the Cauchy problem of the equations of 
radiation hydrodynamics were taken from the stellar evolution calculations for population I 
stars (X = 0.7, Z = 0.02) with initial masses 7OM < M ZA ms < 9OM . This paper continues 
our earlier studies of nonlinear radial oscillations in the massive helium-burning stars (Fadeyev, 
2007; 2008a; 2008b) where the computational methods are described in more detail. 

EVOLUTIONARY MODELS 

In the present study the stellar evolution during the hydrogen core burning was calculated 
with mass loss rates M by Vink et al. (2000, 2001) that are based on the stellar wind models 
with multiple photon scattering and that are in a good agreement with observations of O and 
B stars. Comparison with our previous calculations (Fadeyev, 2007, 2008a, 2008b) shows that 
application of the formula by Nieuwenhuijzen and de Jager (1990) for the hydrogen burning 
phase leads to much higher (by a factor of < 20) mass loss rates and, therefore, to substantially 
smaller stellar mass and luminosity at the core hydrogen exhaustion. For example, evolution 
calculation for the star M ZA ms = 8OM up to the central hydrogen abundance X c = 10~ 5 gives 



the stellar mass and luminosity M = 43.6M Q and L = 1.27 ■ 1O 6 L with mass loss rates by 
Nieuwenhuijzen and de Jager (1990), whereas from calculations with mass loss rates by Vink 
et al. (2000, 2001) the mass and luminosity are M = 50.9M© and L = 1.55 • 1O 6 L . 

The formulae for M by Vink et al. (2000, 2001) become inapplicable after the core hydrogen 
exhaustion and following Vazquez et al., (2007) the mass loss rates at later evolutionary stages 
were calculated according to Nieuwenhuijzen and de Jager (1990). During the helium core 
burning when the effective temperature rises above T e g = 10 4 K the mass loss rates were 
calculated according to Nugis and Lamers (2000). 

Fig. [T] displays the HRD with several LBV-stars with luminositiues and effective temper- 
atures from de Jager (1998). On the same figure are shown the parts of evolutionary tracks 
with initial masses M Z a M s — 7OM and M Z ams — 9OM . The tracks cross the HR-diagram 
from the right to the left and labels at the tracks give the central helium mass fraction Y c . 
Within the displayed parts of the tracks the central temperature T c and the central gas density 
p c change negligibly. For all models T c « 2.15 ■ 10 8 K while the central gas density ranges within 
p c » 250 cm 3 /g for M ZAM s = 9OM and p c « 300 cm 3 /g for M ZAM s = 7OM . 

For M Z ams — 7OM the time needed to cross the effective temperature range 1.5 • 10 4 K < 
T e s < 5 • 10 4 K is ~ 8 • 10 3 yr and the stellar mass decreases from 36M to 31M . For 
^zams = 9OM the evolution time within the same effective temperature range is ~ 5 ■ 10 4 yr 
and the stellar mass decreases from 49M to 36M . Vink and de Koter (2003) evaluated 
the mass of AG Car as M ~ 35M . Bearing in mind existing uncertainties in estimates of 
the luminosities of LBV-stars one may conclude that evolutionary sequences computed in the 
present study are in a good agreement with observations. 

HYDRODYNAMIC MODELS 

Some evolutionary stellar models corresponding to the early helum core burning were used 
as initial conditions in hydrodynamic computations of the self-exciting radial stellar oscillations. 
Main parameters of the hydrostatically equilibrium models are given in the table, where R is 
the equilibrium radius of photosphere. For effective temperatures 1.5 • 10 4 K < T e g- < 3 • 10 4 K 
the mass of the enevelope surrounding the compact core is almost negligible (M env < 10~ 5 M) 
and the mean molecular weight of the stellar material in the envelope does not depend on the 
radius r. Thus, the surface mass fractions of hydrogen X s and helium Y s given in the table are 
the same for all mass zones of hydrodynamical models. 

Outer layers of hydrostatically equilibrium LBV-stars are very close to the boundary of 
dynamical instability due to the large radiation pressure, so that approximation errors of finite- 
difference equations can lead to expansion of outer mass zones with velocity higher than the local 
escape velocity. In such a case the oscillations can be computed only within a too short time 



interval because when the radius of the upper boundary becomes about several dozen times its 
initial value the iteration solution of implicit difference equations do not converge. To overcome 
this obstacle one should diminish approximation errors of the difference equations. After a 
number of test computations it was found that the appropriate solution can be obtained for the 
number of mass zones ranging within 10 3 < iV < 3 • 10 3 . In all hydrodynamic models the size 
of the mass interval decreases to the inner boundary in order to provide enough approximation 
at the core boundary where the gas temperature, pressure and density undergo the sharp rise 
to the center. 

It is assumed that the radius and luminosity at the inner boundary remain constant, that is 
dro/dt = dLo/dt = 0. Determination of the inner radius ro needs the compromise between the 
demands of accuracy and the time step limitations imposed by the Courant stability criterion. 
Choosing the location of the inner boundary and the distribution of mass zones we tried to 
satisfy the condition that the integration time step is At ~ 10 _5 td y n, where tdyn = (R? / GM) 1 / 2 
is the dynamical time scale of the star and G is the Newtonian constant of gravitation. Thus, 
the ratio of the radius of the inner boundary to the equilibrium radius of photosphere is 0.01 < 
Tq/R < 0.05. The gas temperature at the inner boundary ranges within 5 ■ 10 5 K < T < 10 6 K, 
so that therminuclear energy sources can be ignored. 

Hydrodynamic computations with the large number of mass zones are time consuming but 
they allowed us to obtain the solution of the equations of hydrodynamics within time intervals 
as long as 10 4 day and to apply the discrete Fourier transform for determination of the mean 
period II of radial stellar pulsation. It should be noted however that due to irregular dynamical 
behaviour the expansion velocity of the upper boundary sometimes exceeded the local escape 
velocity. In such cases one or a few mass zones excluded from the model and calculations were 
continued with the smaller number of mass zones N. 

All the models of LBV-stars considered in the present study were found to be unstable 
against radial oscillations, the amplitude growth time being comparable with dynamical time 
scale tdyn- The amplitude growth ceases at the amplitude 5r s ~ R, so that radial oscillations of 
LBV-stars are strongly nonlinear. The nonlinear radial oscillations are illustrated in Fig. [2]^a) 
and Fig. E](b) where the temporal dependences of the velocity U s and the radius r s of the upper 
boundary are shown for the model Mzams = 7OM , T e g = 3 • 10 4 K. 

Notwithstanding the large amplitude of the radial displacement the amplitude of light 
changes is < 0.2 mag. Here one should note that the rapid light variations (for the model in 
Fig. [2] with characteristic time < 1 day) are due to the discrete nature of the hydrodynamical 
model and therefore they are not connected with stellar pulsation. 

In contrast to the most classical radially pulsating stars where nonlinear effects become 
important in the outer layers (e.g., W Vir and Mira-type variables) radial oscillations of LBV- 
stars are nonlinear within the entire envelope. This is due to the small pressure gradient and 



approximately constant gas density between the core boundary and the photosphere. Large 
amplitude oscillations in different layers of the pulsating envelope are shown in Fig. [3] where 
plots of the gas flow velocity U are given for several mass zones of the hydrodynamical model. 
For the sake of convenience the plots are arbitrarily shifted along the vertical axis. The lowest 
plot in Fig. [3] corresponds to the layers with the mean radius r < 0.1R. 

The main parameter of the radially pulsating star is the pulsation period II but calculation 
of this quantity for hydrodynamic models of LBV-stars is complicated not only due to irregular 
oscillations. The problem is that the most of the mass is confined inside the compact core and 
the mass of the envelope surrounding the core is M env < 1CT 5 M, whereas its extension is 
about 0.9 of the stellar radius. Because of the low gas density and the small sound speed 
the characteristic motion time of the envelope layers is significantly longer than that of the 
core boundary, so that radial oscillations of LBV-stars are rather the nonlinear running waves 
propagating from the inner core to the stellar surface. Effects of running waves are illustrated 
in Fig. H](a), where three plots of the gas flow velocity U are shown as a function of the 
Lagrangean coordinate measured from the inner boundary. The expansion of inner layers 
occurs approximately at t = t\ and t = t 2 but displacement of the running wave during 
the time interval t?, — t\ is significantly less than the stellar radius. Thus, the running wave 
reaches the upper boundary during several oscillation periods of inner layers. Dependence of 
the oscillation period on the spatial coordinate is clearly seen from velocity plots corresponding 
to different depths inside the envelope (Fig. [3]). 

The bolometric radiation flux from the upper boundary depends on the contribution of 
layers at different depths with different oscillation periods. That is why light variations rather 
weakly correlate with variations of the upper boundary velocity (see Fig. [2]). The pulsation 
period calculated from the discrete Fourier transform of the kinetic energy of the pulsating 
envelope is always longer than that evaluated from the light curve. This is due to the fact 
that the most of the kinetic energy is contributed by outer layers, whereas the light variations 
depend on deeper layers with shorter periods. The typical Fourier spectra of the kinetic energy 
S(E K ) and bolometric light S(AM hol ) are shown in Fig. [5] for the model M Z ams = 8OM , 
T e fj = 2 • 10 4 K. The mean period of variations of the kinetic energy is U(Ek) = 51.2 day while 
the mean light period is II = 16.6 day. 

Last two columns of the table give the mean periods of light changes II and corresponding 
pulsation constants Q. Here one should bear in mind that in contrast to radial pulsations in 
the form of standing waves the pulsation constants of LBV-stars cannot be considered as their 
mechanical characteristics. 

The sharp increase of the luminosity L r at the front of the running wave propagating to 
the upper boundary (see Fig. HI) is the principal cause of the change of radiation flux emerging 
from the upper boundary. Thus, the light variations of LBV-stars are due to dissipation of the 



kinetic energy of running waves rather than due to K-mechanism. Indeed, modulation of the 
radiation flux by K-mechanism can arise in the vicinity of Z-bump (T ~ 2 • 10 5 K) but in LBV- 
stars these layers are at the boundary of the compact core, so that their mass is too small for 
driving the pulsational instability. For some models we carried out hydrodynamical calculations 
with different location of the inner boundary and found that models with temperature at the 
inner boundary T ~ 10 5 K demonstrate instability similar to that obtained for models with 
deeper inner boundary locating below the Z-bump. Therefore, the pulsational instability is 
due to the fact that the adiabatic exponent within the envelope is close to the critical value 
Ti = 4/3. The inefficiency of the /t-mechanism in LBV-stars was considered in the framework 
of the linear theory by Kiriakidis et al. (1993). 

CONCLUSION 

In the present study we computed the hydrodynamical models of radially oscillating LBV- 
stars with effective temperatures 1.5 • 10 4 K < T cff < 3 • 10 4 K and this allows us to conclude 
that microvariability is most probably due to radial oscillations. The K-mechanism is not 
responsible for pulsational instability and oscillations appear is due thermodynamical properties 
of the stellar material with adiabatic exponent close to its critical value T\ = 4/3. Strongly 
nonadiabatic oscillations prevent the development of the dynamical instability due to large 
radiative losses accompanying the motion of the gas. 

Nonlinear radial oscillations of LBV-stars exist in the form of nonlinear running waves 
propagating from the boundary of the compact core to the stellar surface, the amplitude of 
radial velocity variations being as large as several hundred km/s. The amplitude of theoretical 
bolometric light curve is (AM bol < 0.2 mag) seems to be overestimated due to limited spatial 
resolution of hydrodynamical models. 

Mean periods estimated from hydrodynamical models (6 day < IT < 31 day) are in agree- 
ment with observations. For example, observational estimates of microvarion periods are as 
follows: II(P Cyg) « 17.3 day (de Groot et al, 2001), II(B416 M33) « 8.26 day (Shem- 
mer et al, 2000), 20 day < n(v532 M31) < 30 day (Sholukhova et al., 2002), 25.7 day < 
Il(^ 1 Sco) < 32 day (Sterken et al., 1997). The only exception is rj Car but its too long period 
of LT = 58.58 day (Sterken et al., 1996) is rather due to extremely high luminosity (L ~ 5-lO 6 L ) 
of this star. 
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Table 1: Mode 


s of LBV-stars. 


M ZAMS /M Q 


T eS , 10 3 K 
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II, day 


Q, day 
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1.24 • 10 6 
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165 


0.23 
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Figure 1: LBV-stars on the HRD (according to de Jager (1998)) and parts of the evolutionary 
tracks M Z ams = 70 M and M Z ams = 9OM . Attached at the curves are the values of the 
central helium abundance Y c . 
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Figure 2: Hydrodynamical model M ZAM s = 7OM , T eS = 3 • 10 4 K (M = 35.7M , L = 
1.2 • 1O 6 L ). (a) - Gas flow velocity at the upper boundary of the model U s ; (b) - the radius of 
the upper boundary r s in units of the equilibrium radius of photosphere; (c) - bolometric light 
AM bol . 




Figure 3: The gas flow velocity U in some mass zones of the same model as in Fig. El The 
plots are arbitrarily shifted along the vertical axis. Attached at the plots is the Lagrangean 
coordinate 1 — M r /M. 
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Figure 4: (a) - The gas flow velocity U; (b) - the luminosity L r in units of the equilibrium 
luminosity L as a function of the Lagrangean coorinate measured from the inner boundary. 
The hydrodynamics model is the same as in Fig [2] and Fig. [3j In solid, dotted and slashed lines 
are shown the plots for t±, t2 and t%. 
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Figure 5: The hydrodynamical model M ZAM s = 8OM , T cS = 2 • 10 4 K (M — 42.3M , L = 
1.6 • 1O 6 L ). (a) - The Fourier spectrum of the kinetic energy of pulsating envelope S(Ek)', (b) 
- the Fourier spectrum of the bolometric light S(AM^ \). 



